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Abstract Infections by type II feline coronaviruses 
(FCoVs) have been shown to be significantly correlated 
with fatal feline infectious peritonitis (FIP). Despite nearly 
six decades having passed since its first emergence, dif¬ 
ferent studies have shown that type II FCoV represents 
only a small portion of the total FCoV seropositivity in 
cats; hence, there is very limited knowledge of the evolu¬ 
tion of type II FCoV. To elucidate the correlation between 
viral emergence and FIP, a local isolate (NTU156) that was 
derived from a FIP cat was analyzed along with other 
worldwide strains. Containing an in-frame deletion of 442 
nucleotides in open reading frame 3c, the complete genome 
size of NTU156 (28,897 nucleotides) appears to be the 
smallest among the known type II feline coronaviruses. 
Bootscan analysis revealed that NTU156 evolved from 
two crossover events between type I FCoV and canine 
coronavirus, with recombination sites located in the 
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RNA-dependent RNA polymerase and M genes. With an 
exchange of nearly one-third of the genome with other 
members of alphacoronaviruses , the new emerging virus 
could gain new antigenicity, posing a threat to cats that 
either have been infected with a type I virus before or never 
have been infected with FCoV. 

Keywords Type II feline coronavirus • Feline infectious 
peritonitis • Genome analysis • Virus recombination • Virus 
evolution 

Introduction 

Feline coronaviruses (FCoVs) are large, enveloped, posi¬ 
tive-strand RNA viruses with a genome of approximately 
29,200 nucleotides [1-3]. The FCoVs belong to the genus 
Alphacoronavirus , family Coronaviridae , order Nidovi- 
rales. Other members of this subgroup include canine 
coronavirus (CCoV), transmissible gastroenteritis virus 
(TGEV), raccoon dog CoV (RDCoV/GZ43/03), and Chi¬ 
nese ferret badger CoV (CFBCoV/DM95/03) [4]. 

FCoVs are associated with diseases that range from sub- 
clinical and/or mild enteric infections to fatal infectious 
peritonitis [5]. Despite the high prevalence of FCoVs in feline 
populations around the world, only 5-12 % of seropositive 
cats develop feline infectious peritonitis (FIP). FIP is a 
chronic, progressive, immune-mediated disease in domestic 
and nondomestic fields. The typical histopathological finding 
of this disease is systemic perivascular necrotizing pyogra- 
nulomatous inflammation [6]. Two serotypes that differ in 
their growth characteristics in tissue culture and in their 
genetic relationship to CCoV and TGEV have been identified 
[7, 8]. Type II FCoV is significantly correlated with FIP when 
compared to type I viruses [9, 10]. However, unlike the 
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ubiquity of type I FCoV, infection by type II virus encom¬ 
passes only a small percentage of the total number of FCoV- 
seropositive cats in different studies [9-16]. 

Type II FCoV is estimated to have diverged from 
alphacoronavirues in 1953 [4]. Based on partial genomic 
sequence analysis, type II FCoVs were suggested to result 
from a double recombination between type I FCoV and 
CCoV [17]. Like most RNA viruses, Co Vs mutate at a high 
frequency due to the high error rate of RNA polymeriza¬ 
tion. In addition, a unique feature of CoV genetics is the 
high frequency of RNA recombination in the natural evo¬ 
lution of this virus [18]. Recombination among CoVs is an 
attribute of the genus and is thought to contribute to the 
emergence of new pathotypes, such as severe acute respi¬ 
ratory syndrome CoV [19, 20], human CoV NL63 (HCoV 
NL63) [21], HCoV HKU1 [22], and avian infectious 
bronchitis virus (IBV) [23-26]. To gain better evolutionary 
insight into type II FCoVs, we analyzed the complete 
genome of a novel type II FCoV isolate. Taking our data 
together with data from other strains, we discuss the evo¬ 
lution of type II FCoV. 

Materials and methods 

Virus and isolation of viral RNA 

FCoV NTU156 was isolated in 2007 from a kitten with 
naturally occurring FIP by the co-cultivation of pleural 
effusion with feline fcwf-4 cells [27] . After three rounds of 
purification by limited dilution, the virus was propagated 
and titrated. All of the viruses used in this study for the 
sequencing of complete genome came from a stock virus 
passaged six times. NTU156 is relatively fast-growing, 
induces a coronavirus-typical syncytial cytopathic effect 
and is a type II FCoV [10]. 

Reverse transcription and PCR amplification of viral 
sequences 

Eleven microliters of isolated RNA was added to the pre¬ 
mix, consisting of 4 pi of 5 x RT buffer, 2.5 mmol dNTPs 
(GeneTeks Bioscience, Inc., Taipei), 50 pmol random 
primer, 0.2 mol dithiothreitol, and 1 pi of 200 U Moloney 
murine leukemia virus reverse transcriptase (Invitrogen, 
CA, USA) in a 0.6 ml reaction tube. This reaction mixture 
was then briefly centrifuged and incubated at 37 °C for 
60 min, then at 72 °C for 15 min, and finally at 94 °C for 
5 min. A total of 120 primers for PCR were chosen from a 
relatively conserved region of the FCoV genome. Follow¬ 
ing reverse transcription, 1 pi of the RT reaction mixture 
was added to 49 pi of the PCR mixture, which consisted of 
5 pi of 10 x Taq buffer, each primer (10 pmol), dNTP 


(2.5 mmol), 2 U of Taq DNA polymerase (GeneTeks, 
BioScience, Inc., Taipei), and 39 pi of 0.1 % DEPC water. 
An AB1-2720 thermal cycler (Applied Biosystems, USA) 
consisted of 3 min of preheating at 94 °C, followed by 35 
cycles of denaturation at 94 °C for 30 s, annealing at 55 °C 
for 30 s, and extension at 72 °C for 1 min with a final 
extension at 72 °C for 7 min. The viral RNA termini were 
amplified using 3'- and 5'-rapid amplification of cDNA 
ends (Invitrogen, USA). 

Analysis of PCR-amplified products and sequencing 

A total of 10 pi of PCR products from each PCR mixture 
was analyzed using a 1 % agarose gel (Viogene, Taipei) for 
electrophoresis. Amplification products were visualized 
using UV illumination after ethidium bromide staining. 
The nucleotide sequences of the targeted DNA fragments 
were purified (Geneaid Biotech, Taipei) and sequenced in 
both directions using an auto sequencer (ABI 3730XL, 
USA). Full-length genome sequencing of NTU156 was 
performed by single-round PCR with a set of overlapping 
PCR products (average size 750 bp) that encompassed the 
entire genome. 

Sequence analysis 

The complete sequences of NTU156 were then compared 
with other alphacoronaviruses and the results are sum¬ 
marized in Table 1. Multiple alignments of nucleic acid 
sequences were performed by the Clustal W method using 
the Meg Align program (DNASTAR Inc., WI, USA). 
Phylogenetic analyses were conducted using MEGA, ver¬ 
sion 4.0. Similarity graphs were prepared with SimPlot 2.5 
software [28]. Potential recombination sites were identified 
using the Recombination Detection Program (RDP) [29]. 


Results 

Genomic sequence of FCoV NTU156 

The full genomic RNA sequence of FCoV NTU156 com¬ 
prises 28,897 nucleotides (nts), excluding the 3' poly- 
adenylated nts. Sequence analysis revealed that NTU156 
contains conserved open reading frames with an overall 
genome organization similar to known FCoVs (Table 2). 
The overall nucleotide composition is as follows: A, 
29.3 %; C, 17.3 %; G, 21.0 %; and T, 32.4 %. The G+C 
content is 38.3 %. NTU156 possesses the putative tran¬ 
scription regulatory sequence (TRS) motif, 5 / -CUAAAC-3 / , 
at the 3' end of the leader sequence and preceding each ORF 
(Table 2). 
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Table 1 Isolation dates and 

accession 

numbers of NTU156 and other alphacoronavirses 





Isolation 




Virus 

Host 


Year 

Country 

References 

GenBank accession no. 

FCoV NTU156 a 

Cat 


2007 

Taiwan 

Fin [10] 

GQ152141 

FCoV NTU2 

Cat 


2003 

Taiwan 

Unpublished 

DQ160294 

FCoV 79-1146 a 

Cat 


1979 

USA 

McKeirnan (1981) 

DQ010921 

FCoV Black a 

Cat 


1980 

USA 

Black (1980) 

EU186072 

FCoV UU2 a 

Cat 


1993 

USA 

Unpublished 

FJ938060 

FCoV DF-2 a 

Cat 


1981 

USA 

Everman (1981) 

JQ408981 

FCoV ClJe a 

Cat 


2007 

UK 

Dye [2] 

DQ848678 

FCoV 79-1683 a 

Cat 


1979 

USA 

McKeirnan (1981) 

JN634064 

FCoV KU-2 

Cat 


1991 

Japan 

Hohdatsu (1991) 

D32044 

FCoV Yayoi 

Cat 


1981 

Japan 

Hayashi (1981) 

AB695067 

FCoV C3663 

Cat 


1994 

Japan 

Mochizuki (1997) 

AB535528 

CCoV NTU336 a 

Dog 


2008 

Taiwan 

Unpublished 

GQ477367 

CCoV 1-7 l a 

Dog 


1971 

Germany 

Binn (1974) 

JQ404409 

CCoV TN-449 a 

Dog 


1980s 

USA 

Unpublished 

JQ404410 

CCoV 174/06 

Dog 


2006 

Italy 

Decaro (2009) 

EU856362 

TGEV Purdue a 

Pig 


1946 

USA 

Doyle (1946) 

DQ811789 

TGEV M6 a 

Pig 


1965 

USA 

Bohl (1965) 

DQ811785 

TGEV TS a 

Pig 


2004 

China 

Cheng (2004) 

DQ201447 

HCoV 229E a 

Human 


1963 


Hamre (1966) 

AF304460 

a Full genome available 







Table 2 Coding potentials and putative transcription regulatory sequences (TRS) of FCoV NTU156 


ORFs Start-end (nucleotide 

No. of nucleotides 

No. of amino acids 

Putative TRS 


position) 












Nucleotide position 

TRS sequences 

ORF la 313-2,417 


12,105 


4,035 

94 

GAACUAAACGAA210 AUG 

ORF lab 313-20,417 


20,105 


6,701 



(shift at 12,381) 






S 20,414-24,775 


4,362 


1,453 

20,382 

UUACUAAACUUU23AUG 

ORF 3a 24,840-25,076 


237 


78 

24,832 

GAACUAAACUUAUG 

ORF 3b 25,021-25,236 


216 


71 



ORF 3c 25233-25,433 


201 


66 



E 25,512-25,760 


249 


82 

25,469 

GUU CU A A ACG A A3 7 AU G 

M 25,771-26,556 


786 


261 

25,762 

GAACUAAACAAAAUG 

N 26,569-27,687 


1,119 


372 

26,557 

UAACUAAACUUCUAAAUG 

ORF 7a 27,692-27,997 


306 


101 

27,684 

GAACUAAACGCAUG 

ORF 7b 28,002-28,622 


621 


206 




5', y Untranslated regions (UTRs) and open reading 
frames (ORFs) 

The 5'-UTR of FCoV NTU156 comprises 312 nts—one or 
two nts more than other known FCoVs (79-1146: 311 nts, 
black: 311 nts, ClJe: 310 nts)—whereas the 3 / -UTR 
of our virus comprises 275 nts, which is identical to other 
FCoVs. The 5' two-thirds of the genome contain the la 


(nt 313-12,417) and lab (nt 313-20,417) genes that encode 
the nonstructural polyproteins. A typical coronavirus “slip 
site,” 5 / -UUUAAAC-3 / (nt 12,381-12,387) is located within 
this gene. These genes are followed by genes encoding the 
four structural proteins: spike (nt 20,414-24,775), envelope 
(nt 25,512-25,760), membrane (nt 25,771-26,556), and 
nucleocapsid (26,569-27,687) (Table 2). The two strings of 
accessory genes identified in all of the known FCoVs, 
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i.e., ORF 3ab and ORF 7ab, were found in NTU156 as well 
(Table 2). However, an in-frame deletion of 442 nucleotides 
in ORF 3c was identified, which resulted in a relatively short 
gene comprising only 201 nts. 

Genetic comparison and phylogenetic analysis 
with other alphacoronaviruses 

The overall sequence comparison revealed that NTU156 
was more closely related to known subgroup la CoV but 
not lb within alphacoronaviruses (Fig. 1). Nucleotide 
sequences similarity graphs of NTU156 with known type I 
FCoVs, CCoVs, and TGEVs were created by the Simplot 
software. The results showed that NTU156 was more clo¬ 
sely related to type I FCoVs from the 5' end of the genome 
to position 15,000 and from position 27,500 to 3 / -UTR 
(Fig. 2). Genes located at the 5' end (nspl-11) and 3' end 
(the N gene through ORF7) of NTU156 show consistently 
high similarity to type I FCoVs, whereas from nspl3 
through the E sequence, the similarity to canine and por¬ 
cine Co Vs varies dramatically (Fig. 3). These data indicate 
that NTU156 might have arisen from recombination events 
between different strains of Co Vs from species other than 
cats. Two possible recombination sites, at approximate 
positions 14,300 and 27,300, corresponding to the RNA- 
dependent RNA polymerase (RdRp) and the M gene, 
respectively (Fig. 4), were further analyzed. 

Phylogenetic trees using the nucleotide sequence of 
genes for putative proteins and polypeptides of alphacor¬ 
onaviruses were further constructed. At nsp 1 through nsp 
11 (Supplementary Fig. la) and from N gene (Supple¬ 
mentary Fig. lb) to the ORF 7ab gene (Supplementary 
Fig. lc), NTU156 was more closely related to type I 
FCoVs. At the nsp 12 (RdRp) and the M gene, NTU156 
was not clustered with any known alphacoronavirues 
(Fig. 3a, f). From the nspl3 through the E gene, NTU156 



Fig. 2 Nucleotide sequence similarity of the complete genome of 
FCoV NTU156 with other alphacoronaviruses. A similarity plot was 
constructed to identify the sequence homology between type I FCoVs 
black, ClJe, and UU2 (gray); CCoV NTU336 (red); and TGEV 
purdue, M6, and TS (blue). Red arrows represent putative recombi¬ 
nation regions. A similarity of 1.0 indicates regions that share 100 % 
nucleotide identity. The similarity calculation was performed using 
the following parameters: a window size of 1,000 bp and a step size 
of 200 bp for the full-length sequence 

was clustered with CCoV (Fig. 3b-e). Taken together, 
these data indicated that NTU156 might have evolved from 
two recombination events with CCoV, with the sites of 
recombination located in the RdRp and M genes. 

Discussion 

A unique feature of CoV genetics is the high frequency of 
RNA recombination both in vivo and in vitro [18]. Here, an 
interspecies recombination between feline and canine CoV 
was identified in a viral strain NTU156, which was isolated 
from the pleural effusion of a FIP cat. This is the first time 
that evidence for natural recombination has been docu¬ 
mented through the complete genome sequence analysis of 
type II FCoV. In 1998, Herrewegh et al., based on partial 


Fig. 1 Phylogenetic 
relationships constructed using 
the complete genome sequences 
of FCoV NTU156 and other 
alphacoronaviruses . Analysis 
was performed using MEGA 4 
software and neighbor-joining 
methods based on 1,000 
replicates. Bootstrap support 
values greater than 90 are 
shown 
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Fig. 3 Phylogenetic relationships constructed using nsp 12 (a), 
nsp 13-16 (b), S gene (c), ORF3 gene (d), E gene (e), and M gene 
(f) sequences of FCoV NTU156 and other alphacoronaviruses. 


Analysis was performed using MEGA 4 software and neighbor¬ 
joining methods based on 1,000 replicates. Bootstrap support values 
greater than 90 are shown 


sequence analysis, first determined that type II FCoVs 
79-1146 and 79-1683 originated from a homologous RNA 
recombination event between type I FCoV and CCoV [17]. 
The complete sequence of strain 79-1146 was later pub¬ 
lished in 2005 [1]. When comparing strain 79-1146, 
79-1683, DF-2, and our NTU156 strain, the only four type 
II FCoVs that have had their full-length genomes 
sequenced to date, a common phenomenon was found; both 
viruses arose from a double recombination event between 
type I FCoV and coronaviruses from other species. 

NTU156 appears to have evolved from a recombination 
between type I feline and canine CoV; however, when we 
aligned the genes in which recombination took place with 


other type II FCoVs, genome crossovers with other 
alphacoronavirus were noted. When the sequences around 
the putative recombination sites were examined, i.e., one 
located in the 5' region (strain 79-1146) and two in the 3' 
region (strain 79-1146, and 79-1683), porcine coronavirus 
(TGEV) was also found to contribute to the evolution of 
type II FCoV (in addition to CCoV) (Figs. 3b, 5). This 
finding is not surprising because the receptor for type II 
FCoV, feline aminopeptidase N, has been found to serve as 
a receptor for several CoVs, including canine, porcine, and 
human CoVs [30]. Therefore, the interspecies recombina¬ 
tion of type I FCoV with any of the above viruses might 
occur in nature. The recombination of type I FCoV in cats 
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Fig. 4 Bootscan analysis of the complete genome of FCoV NTU156 
with other alphacoronaviruses (type I FCoV UU2: pink ; CCoV 
NTU336: yellow ; TGEV purdue: green). Analysis was performed 
using RDP software and neighbor-joining methods based on 1,000 

living in the same household with dogs or living close to 
pig farms could give rise to type II FCoV. 

Based on the analysis of the four full genomes of type II 
FCoV available at present (NTU156, 79-1146, 79-1683, 
and DF-2), type II FCoVs appear to retain type I FCoV 
sequences in their 5' and 3' ends. We asked whether the 
genes located in these regions are indispensable for FCoV 
replication in cats. To answer that question, the amino acid 
sequence of genes retained at both ends of FCoVs, CCoVs, 
and TGEVs, i.e., nsp 1 through nsp 11 and the N and ORF7 
genes, were further aligned and compared (Table 3). In 
contrast to the greater than 90 % of amino acid homology 
between different strains of FCoVs, the nsp 2, 3, and 6, as 
well as the N gene and the ORF7a gene of FCoVs, when 



Fig. 5 Nucleotide sequence similarity of the complete genome of 
FCoV 79-1146 with other alphacoronaviruses. A similarity plot was 
constructed to identify the sequence homology between type I FCoVs 
black, ClJe, and UU2 (gray); CCoV NTU336 (red); and TGEV 
purdue, M6, and TS (blue). Red arrows represent putative recombi¬ 
nation regions. A similarity of 1.0 indicates regions that share 100 % 
nucleotide identity. The similarity calculation was performed using 
the following parameters: a window size of 1,000 bp and a step size 
of 200 bp for full-length sequences 


replicates. Putative recombination regions are indicated with arrow. 
The boo strap support calculation was undertaken using a window size 
of 1,000 bp and a step size of 200 bp. Dashed line indicates the 95 % 
bootstrap support threshold values 

compared with CCoVs or TGEVs, exhibited similarity of 
less than 80 %. This finding indicates that these gene 
products might possess irreplaceable functions for FCoV 
replication. This might explain why type II FCoVs found in 
nature harbor genomes that evolved from a double 
recombination. 

Although the prevalence of type II FCoV is consistently 
lower (2-11 %) than type I virus around the world 
(88-98 %) [9-16], our previous study indicates that 
infection of type II FCoV correlates significantly with FIP 


Table 3 Percent similarity with respect to the amino acid sequences 
of FCoV NTU156 and other alphacoronaviruses 


Gene 

FCoV 



CCoV 


TGEV 



NTU2 

Black 

ClJe 

NTU336 

174/06 

Purdue 

TS 

ORF1 

nspl 

a 

93.6 

92.7 

91.8 

— 

92.7 

91.8 

nsp2 

— 

93.6 

91.7 

80.2 

— 

79.1 

78.7 

nsp3 

— 

91.5 

84.9 

79.2 

— 

78.8 

78.3 

nsp4 

— 

96.3 

93.5 

88.8 

— 

88.2 

87.1 

nsp5 

— 

97.4 

96.7 

93.4 

93.4 

93.0 

92.7 

nsp6 

— 

93.9 

91.5 

77.9 

76.5 

77.9 

77.2 

nsp7 

— 

97.6 

97.6 

97.6 

96.4 

96.4 

96.4 

nsp8 

— 

97.9 

96.4 

92.8 

90.8 

92.8 

92.3 

nsp9 

— 

95.5 

96.4 

86.5 

91.0 

89.2 

89.2 

nsp 10 

— 

98.5 

97.8 

95.6 

95.6 

95.6 

94.8 

nspll 

— 

94.7 

100 

84.2 

84.2 

84.2 

84.2 

N 

93.3 

92.0 

91.7 

77.5 

77.2 

76.4 

76.9 

ORF7 

7a 

99.0 

97.1 

96.1 

81.4 

81.4 

78.5 

74.7 

7b 

88.9 

86.9 

86.7 

59.9 

61.8 

— 

— 


a Not available 
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when compared to type I [10]. As shown in the present 
study, type II FCoV arises by exchanging a large genome 
fragment (approximately 12 kb) of type I FCoV with other 
members of alphacoronaviruses. The genes exchanged 
through this double recombination include nsp 13-16, 
structure protein S (spike), and accessory protein 3abc. The 
nsp 13-16 proteins are replication proteins with functions 
such as helicase activity (nsp 13), nucleoside triphospha¬ 
tase activity (nsp 13), RNA S'-triphosphatase activity (nsp 
13), 3'-5' exoribonuclease activity (nsp 14), RNA cap 
formation (nsp 14 and nsp 16), and endonuclease activity 
(nsp 15) [31]. It has ben reported that the function of the 3c 
protein might be crucial for viral replication in the gut but 
is dispensable for systemic FCoV replication [32], How¬ 
ever, S proteins play a crucial role in receptor binding and 
eliciting protective immunity [33]. Through the replace¬ 
ment of nearly one-third of the genome, the new virus 
might gain new antigenicity, posing a threat to cats that 
either have been infected with a type I virus before or never 
have been infected with FCoV. 
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